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We study the evolution of the superconducting coherence length from weak to strong coupling, 
both within a s-wave and a d-wave lattice model. We show that the identification of with the 
Cooper-pair size ^pair in the weak-coupling regime is meaningful only for a fully-gapped (e.g., s- 
wave) superconductor. Instead in the d-wave superconductor, where ^pair diverges, we show that 
^0 is properly defined as the characteristic length scale for the correlation function of the modulus 
of the superconducting order parameter. The strong-coupling regime is quite intriguing, since the 
interplay between particle-particle and particle-hole channel is no more negligible. In the case of 
s-wave pairing, which allows for an analytical treatment, we show that is of order of the lattice 
spacing at finite densities. In the diluted regime diverges, recovering the behavior of the coherence 
length of a weakly interacting effective bosonic system. Similar results are expected to hold for d- 
wave superconductors. 



PACS numbers: 74.20.-z,74.20.Fg,71.10.Fd 

I. INTRODUCTION 

In the last years the analysis of the phase-fluctuation 
contribution to the low-temperature properties of high- 
Tc superconducting cuprates stimulated the interest in 
the microscopic derivation of the phase-only action for 
d-wave superconductors The symmetry of the or- 

der parameter has been shown to play a crucial role in 
determining the quantum-to-classical crossover for phase 
fluctuations via enhanced dissipative effects The 
effective hydrodynamic phase-only action is usually cut 
off at large momenta. A reasonable spatial bound for the 
phase-only action is the coherence length which sets 
the distance above which the fluctuations of the modulus 
of the order parameter |A| become uncorrelated. Indeed, 
at distances greater than the relevant degrees of free- 
dom are the phase fluctuations only. Within this context, 
a proper definition of must be related to the charac- 
teristic length scale of the spatial decay of the correlation 
function for |A|, whatever is the symmetry of the order 
parameter. 

Another relevant length scale for superconductors is 
the characteristic size of the Cooper pair 

_ / /dr|^(r)Pr2 _ / / dk [c?,0(k)]^ 

where V'(r) is the Cooper-pair wave function, 9k is the 
gradient operator in k-space, and (/)(k) = Ao7k/i?k is 
the Fourier transform of V'(r) (|]- Here Aq is the am- 
plitude of the superconducting gap, 7k is the factor 
which controls the symmetry of the BCS wave function, 
i?k = v^k+^o7k' is the normal-state dispersion. 

Throughout this paper we assume the lattice spacing a 
as the unit length. 

It is known that at weak coupling for a s-wave super- 
conductor (i.e., 7k = 1) both and £,pair are, within a 



numerical factor, of the same order vp/Ao Q. Here vp is 
a properly defined average value of the quasiparticle ve- 
locity at the Fermi surface (see Sec. II), which reduces to 
kp/m in a continuum model. By analogy, one would ex- 
pect the same result to hold for d-wave superconductors 
(7k = 7k = coskx — cos ky). However, when evaluating 
Eq. (|l|), one finds that the mean size of the Cooper pair 
is infinite. This fact, which has been often overlooked 
(see, however, Ref. 0), can be easily seen by evaluating 
the dominant contribution to Eq. (|l]), which comes from 
the nodal quasiparticles. As we show in Appendix A, the 
nodal quasiparticles give a logarithmically divergent con- 
tribution to Eq. (|l|). The question then arises whether 
^0 is well defined in the case of d-wave pairing, and of the 
expected order of magnitude in the weak-coupling limit. 

The strong-coupling limit is even more intriguing, as 
one would expect that the magnitudes of the two lengths 
are widely separated. For instance, in the case of s- 
wave pairing, while £,pair should decrease as the pairing 
strength increases, should resemble the behavior of 
the natural length scale for the spatial variation of the 
condensate wave function of an effective "bosonic" sys- 
tem [HJoj. Thus ^0 should be controlled by the inverse 
of the residual interaction between the bosonic Cooper 
pairs, and could eventually increase as the residual boson- 
boson interaction decreases and by approaching the di- 
luted limit. Such a separation between the two length 
scales has indeed been found in Ref. within the con- 
tinuum s-wave model a la Nozieres and Schmitt-Rink 

In this paper we address the issue of the system- 
atic analysis of the behavior of the coherence length 
^0 defined through the spatial decay of the static cor- 
relation function for |A|, A"A(r), which is obtained by 
Fourier-transforming the correlation function in momen- 
tum space A'A(q). We specifically consider the two- 
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dimensional negative-C/ Hubbard model for s-wave pair- 
ing and its simplest generalization for d-wave pairing. 

We show that, whatever is the symmetry of the order 
parameter, in the weak-coupling limit is finite and of 
the expected order vp/ Aq. 

In the strong-coupling limit the modulus of the order 
parameter |A| and the density of particles p are coupled. 
Specifically, in the case of s-wave pairing, they experi- 
ence the same fluctuations: in particular at low density 
the two fields become proportional As a con- 

sequence, the density fluctuations contribute to X^{q) 
via the gap-density coupling. We derive X/\{q) in s- 
wave superconductors by including density fluctuations 
(whose contribution to is negligible at weak coupling) 
and investigate the length scale of the spatial decay of 
its Fourier transform X/^ (r) by varying the density p be- 
tween and 1 (the range 1 < p < 2 being recovered 
by particle- hole symmetry). In the diluted limit {p « 0) 
Xa{^) decays exponentially with a length scale diverg- 
ing as 1/ ^/p, as it is expected for a weakly-coupled diluted 
Bose liquid |l2|J|,^. At higher densities, and for strong 
coupling, A'A(r) oscillates with the periodicity a of the 
lattice. This is due to the fact that XA{q) is dominated 
by momenta q « Q = (tt, tt), in analogy with the density 
mode, which is massless for q = Q at half flUing . By 
means of numerical calculations we show that the length 
above which ArA(r) is strongly suppressed with respect 
to its value at r = is of the order of the lattice constant 
for all densities away from p ~ 0, even though it exhibits 
a long-living exponential tail governed by an increasing 
characteristic length scale approaching half filling. 

The case of d-wave pairing in the strong-coupling limit 
would be much more difficult to address, as the coupling 
of the gap fluctuations to the particle-hole channel re- 
flects in a Hartree-Fock-like correction to the bare band 
dispersion, making the analytical treatment not viable. 
Despite these complications, however, the physics is not 
expected to be different from the s-wave case, which al- 
lows for a much simpler and transparent analytical treat- 
ment. For this reason, we devote most of our strong- 
coupling analysis to the s-wave case. 

Finally, we comment on the consequences of such an 
analysis on the problem of a proper definition of the 
momentum cutoff for phase fluctuations in the effective 
phase-only model. 

The plan of the paper is the following. In Sec. II 
we discuss the weak-coupling regime for both s- and d- 
wave superconductors. In Sec. Ill we devote our analysis 
to the strong-coupling regime for s-wave superconduc- 
tors. In particular. Sec. Ill A addresses the low-density 
regime. Sec. Ill B deals with the intermediate- and high- 
density regime, while in Sec. Ill C we discuss the proper- 
ties of the density-density correlation function. Conclud- 
ing remarks are found in Sec. IV. Detailed calculations of 
the coherence length and of the Cooper-pair size for both 
s- and d-wave superconductors at weak coupling are re- 



ported in Appendix A. In Appendix B we discuss how the 
strong-coupling results discussed in Sec. Ill are affected 
by the inclusion of a next-to-nearest-neighbor hopping 
term. 



II. THE WEAK-COUPLING REGIME 

We start with the BCS action on a two-dimensional 
lattice, at a temperature T — P^^- 

^ = ly^^Y. ^k'x (r) (dr + a)ck. (r) + Hj (r) | , 

k,k',q 

Here Ng is the number of lattice sites, ?7 > is the pair- 
ing interaction strength, — — 2i(cos + cos ky) ^ p is 
the band dispersion associated with a nearest-neighbor 
hopping i, p is the chemical potential, and the factor 7k 
controls the symmetry of the gap. To derive the cor- 
relation function for |A|, we first perform the standard 
Hubbard- Stratonovich decoupling of Hj and then make 
explicit the dependence on the phase 6 and on the mod- 
ulus |A| of the complex order parameter A = |A|e*^ by 
means of the gauge transformation c„ c^e^^^^ p|,p|,p^ . 
As a consequence, the action depends on |A| through the 
Hubbard-Stratonovich Gaussian term and the interaction 
of the field with the fermions, which arises from Hi, i.e. 



5(|A|)=^i|A|,|A 



T ( |AL._, 







where q = (q, ilm), k = (k, Wn), and ojn are 

the bosonic and fermionic Matsubara frequencies re- 
spectively, and we introduced the Nambu spinor 'J'^ = 
(c^|, c_fc,x). We assume, as usual, small fluctuations 
of |A|g around its saddle-point value Aq, i.e., |A|g = 
Aq + (5|A|q. After integrating out the fermionic degrees 
of freedom around the BCS superconducting saddle-point 
solution, we expand the resulting effective action Sef / for 
(5|A|q up to the second order, obtaining 



5e//(<S|A|)=^ 

q 

T 



1 1 , ^' 



^|A|,5|A|_„ 



(3) 



D{q, Qr, 



E^^Tr 



^0 (k + -,u;„ 



X Tl 



.o(^-f. 



where Qo is the Nambu Green function evaluated at the 
BCS level, ri is the Pauli matrix and the trace is taken 
in the Nambu space. 
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We define 



where the "mass term" is given by 



^A(q,a„) = 2<(5|A|,^|A|_, >, 



(4) 



i.e., the inverse of the coefheient of the Gaussian term in 
Eq. (H). Since we are interested in the spatial variation 
of the static correlation function for (5|A|q at zero tem- 
perature, we evaluate the T = limit of XA(q, i^m — 
0) = XA(q), 



^A(q) 
D 



1 1 



U 2 



^(q) 



Ns'^ E+ + E- 



(5) 

l + ^±^-^^l, (6) 



E+E- 



where, with a standard notation. 



Ao7k 



A?,, and A±, E± are calculated at mo- 



t2 

^k ^ "k' 

menta k ± q/2 respectively. 
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FIG. 1. X^^(q) at intermediate coupling I] jt = 1 in the 
case of d-wave pairing for some values of the density p. The 
wave vector q varies along the diagonal of the Brillouin zone, 
and an analogous behavior is observed by varying q along 
different directions. As it is expected in the weak-coupling 
regime, X^^(q) exhibits a minimum at q = [which corre- 
sponds to a maximum for the correlation function XA(q) in 
Eq. §]. 

The length scale which controls the long-distance be- 
havior of Xi^[v) may be extracted from the dominant 
part of its Fourier transform XA(q). It is generally ex- 
pected that the main contribution to Xi^(c\) comes from 
the region q w 0. This is confirmed by our numerical 
calculations. In Fig. |^ we show the q dependence of 
^A^(q) various densities p, for U/t = 1, in the case of 
d-wave pairing. Similar results are found in the s-wave 
case. Expanding Eq. (||) for small momenta, we get 



^A(q) 



1 



mP- -\- cq^ 



(7) 



m 



-E 



2El 



(8) 



and 



= E {^CkAk(2A2 - UDich.S,^) ■ (akAk) 



^[Ck(C^-2A^)a2a + 3Ak^2g^Ak] 



(9) 



When the relevant fluctuations are those near q = 0, 
^A(q) is maximum at q = 0, and c > 0. In such a case 
the approximate expression (|^) can be used to determine 
the long-distance behavior of Xa{j^)- Indeed the Fourier 
transform of Eq. (Q) has an exponential decay in real 
space with a characteristic length scale 



(10) 



?o - \/ — ^ 



As we show in Appendix A, differently from ^pai,- given 
by Eq. (|l|), such a length scale, evaluated from Eqs. (H) 
and (^, is finite also in the presence of nodal quasiparti- 
cles, making it a natural candidate for the characteristic 
length of c?-wave superconductors. In the weak-coupling 
regime {U/t <C 1), ^o, Eq. ([To|), can be expressed both 
in the s-wave and in the d-wave case by means of the 
average values of the Fermi velocity (vp) and of the gap 
{Ap) on the Fermi surface. 



Co 



VF 



12Af 



(11) 



Here we define the average value of a function on the 
Fermi surface as 



{hl)F 



Ekfe^'^(Ck) 

Ek^(ek) 



(12) 



In the s-wave case (7k — 1) Ap — Aq, and ^0 coin- 
cides within a numerical factor of order one with the 
Cooper-pair size (see Appendix A). In the d-wave case 
Ap — AqTp, so that the effective gap A^ which appears 
in Eq. (O) is smaller than the maximum value of the 
gap at the Fermi surface. In both cases, however, in the 
weak-coupling regime, the correlation function for (5|A|q 
exhibits an exponential decay over a length scale of order 
vp/ Aq, while the different gap symmetries only introduce 
a numerical factor jp . It is then natural to assume such a 
length scale as the spatial cutoff for phase fluctuations in 
the effective phase-only action, both in s- and in d-wave 
superconductors. 
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III. THE STRONG-COUPLING REGIME 

The extension of the above resuhs to the strong- 
couphng regime U /t ^ 1 is quite intriguing. Indeed, as 
the pairing increases the coefficient c of Eq. (^) decreases, 
and becomes negative, making the definition ( p^ ) mean- 
ingless. On the other hand, it is commonly expected that 
in the diluted regime the fermionic system maps into a 
(weakly interacting) bosonic system, where the Cooper 
pairs act as boson particles, with a weak residual repul- 
sion between them It is well known that in a 
Bogoljubov liquid of weakly interacting bosons, the co- 
herence length ffjos which controls the correlations of the 
superfluid order parameter diverges as the density de- 
creases. According to Ref. M 



1 



(13) 



where g is the local repulsion, ruB and pB are the mass 
and the particle density of bosons respectively. One 
would like to find a similar behavior for within the 
Hubbard model at strong coupling and in the low-density 
limit. As we anticipated in Sec. I, in the following we 
devote our attention to the case of s-wave pairing, which 
allows for an analytical treatment. However, one expects 
similar results also in the case of d-wave pairing. 

At strong coupling the definition ( [lO| ) must be general- 
ized to include two effects: (i) the fluctuations of |A| are 
strictly tied to the fluctuations of p, in such a way that 
both contribute to the correlation function for |A|; when 
this interplay is considered, we find that (ii) the domi- 
nant contribution to Xa{(1) arises from q ^ (except 
for the small-density regime), and the small- momentum 
approximation (|^) is no longer appropriate. 

We first address the point (i). In the s-wave Hubbard 
model the interaction Hj can be decoupled both in the 
particle-particle and in the particle-hole channel. When 
the density fluctuations are taken into account on the 
same footing as the Cooper-pair fluctuations, the effec- 
tive action (g) gets modified, and reads 



[1=1,2 



(14) 



where <I>|^-^ = bpq,^ 



kXpiq), Xr2^{<l) = X^i^^q) = XpA{q), and X^^^{q) = 
1/U — ^D{q). Below, we only need the expressions in 
the static limit. 



Xp(q) = 



-E 



1 



1 



^0 



E+E^ 



E+ + E- E+E- 



(15) 
(16) 



in Eq. (|T^, and recover the action for 5\A\ only. The 
correlation function now reads 



v - lXp{<l) 



(17) 



It is worth noting that integrating out the density at 
the Gaussian level corresponds to performing the RPA 
resummation for the correlation function in the 
particle-hole channel. At weak coupling the bubble xpA 
which couples the two channels is negligible, and the re- 
sult (|^) is recovered. In order to estimate Eq. (|l^) at 
strong coupling, we evaluate the q-dependent leading or- 
der in t/[/ of the bubbles (|), (|l|), and i^: 



D{ci) ~ D{0) ~ CAw(q) 

-i-E^-^A^q), 



XpA(q) ^ Xpa(O) - CpAw(q) 



E 



CpAw(q), 



Xp(q) 



Xp(0) - Cpw(q) 



(18) 



(19) 



(20) 



where the function i/;(q) = 4(sin^ ^ -|-sin^ ^) reduces to 
q'^ at small momenta, while respecting the lattice peri- 
odicity at higher momenta. In Eqs. (p^)-(20) the q = 
values are of order l/?7, while the coefficients ca, CpA and 
Cp are of order jU^ . All these coefficients are functions 
of the gap amplitude Aq and of the chemical potential /i, 
whose dependence on the density can be determined by 
using the self-consistent saddle-point equations 



1 = 



(5 = 



U ^ 1 



N. 



k 



where we introduced the "doping" 5 = 1 — p to simplify 
the notation [|l^. The above equations can be explicitly 
solved at strong coupling giving 



Ao 



V 



(l-2a2)^l_^2+0(^4^ 



^ = -^[<5(l + 4a2) + o(«4)] 



To take into account the effect of density fluctuations 
on Xa we integrate out the density-fluctuation field 



where we introduced the small parameter a = 2t/U . 
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As a consequence, for a S we can put 

1 



1 



1 1, 



1 - 



icpw(q) 



U - 3Xp(0) 
Substituting Eq. (|2|) into Eq. (|4|) we have 

^A(q) = ^ 



M2 + Cw(q) ^ 



where the q = value = ^ (0) is given by 



(25) 



(26) 



FIG. 2. Doping dependence of the coefficients ca, CpA and 
Cp given by Eqs. (ji^), (|2|) and (^ respectively. 

As anticipated above, we shall investigate the region 
< p < 1, i.e. < S < 1, while the range 1 < p < 2 can 
be recovered by particle- hole symmetry {S ^S). 

With long but straightforward calculations it can be 
seen that at leading order in a 



CA 



+ 1), 



Cp = ^(l-<52)(552-l), 



(21) 
(22) 
(23) 



whose behavior is reported in Fig. |^. 

By means of Eqs. (|ll)-(|2^) and Eqs. (|ll)-(il) we ob- 
tain the q-dependcnt strong-coupling expression of X/\ 
at leading order in a. 



1 1 

U ^2 



(24) 



[DiO) - CA«;(q)] 

X^a(O) 
[Xp(0) - Cpw{q)] 

2XpA(0)CpAw(q) 

^-i[Xp(0)-CpU;(q)]' 
where we neglect the subleading c^^ term. 

A. Low-density regime 



In the limit of low density ((5 w 1) it is possible to 
rewrite Eq. (U) in a simpler form. Indeed, at the lead- 
ing order in a the denominator of the third and fourth 
terms in Eq. is given by 



jj-^ [Xp(0) - Cpu;(q)] ~ ^ 



ia - Yw(q) 



U 

6\1 - 36a' + GOa'S') 



[l - 6a' + m~a 



S'+4a'{l-f6' + fS^) 



'a' 
Aa' 1 - 5' 



and 



1-5(5^ + 5(5 

5'{h5'-l){l-5'f 
((52 + Aa'f 



2 , , 25'{l-5')[55'-3) 



(52 + 4a2 
2a2 25' - 1 



(27) 



(28) 



Therefore, C > in the diluted limit ((5 w 1), and 
A'A(q) is dominated by the small-momentum region, 
where w(q) ~ q^ and the coherence length is given by 
the generalization of Eq. ([lO|), with c and m' substituted 
by C and M' respectively, 



c 

Af2 



I 25' -I 

8(1-^2) 



4Vp 



(29) 



The result (^9|) for the coherence length shows, at low 
particle density, the same divergence of the bosonic co- 
herence length ^bos given by Eq. ( [l3|). Notice that a 
more strict comparison between Eq. (|29|) and Eq. ( |T^ ) 
requires a dependence of on the effective mass of the 
electron pair mp and on the pair-pair residual repulsion 
(7p in the bosonic limit, of the form ~ l/(gpmpp). 
This is indeed the case, since in the bosonic limit of the 
fermionic model pB = p/2, mp ~ U/8t', and the residual 
repulsion of the bosonic model corresponds to the inverse 
of the compressibility of the fermionic system gp ~ x~^- 
Following the analysis of Ref . , it can be seen that at 
strong coupling the compressibility of the Hubbard model 
is X — U/ 8t' jl^ , leading to the explicit cancellation of t 
and [/ in Eq. (13), in agreement with Eq. (p9|). This sup- 
ports the reliability of the mapping of the diluted large-C/ 
Hubbard model into an effective boson model p^ , p^ . 
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B. High-density regime 



As S decreases, according to Eq. (pSD, the coefficient 
C decreases and becomes negative for S < l/\/2- This 
change of sign is not due to a failure of our approxima- 
tion, since at strong coupUng a <^ 1, so that at S = 1/^/2 
the simphfied expression (|6|) is still valid. The vanish- 
ing of the leading coefficient in the small-q expansion 
is instead a signal of the fact that the maximum of the 
correlation function Xa (q) is moving to a finite q. Since 
the q-dependence of Xa (q) comes entirely from the func- 
tion w(q), the only candidate alternative to q = (0,0) is 
q = Q = (7r,7r). Indeed, evaluating = X'^^{Q) ac- 
cording to Eq. (M) we find that 



Mi 



4a^ 
U 



(30) 



at all doping, so that < as far as (5 < l/\/2, as 
depicted in Fig. ^ In such a case, the dominant S\A\q 
modes are located near the wave vector Q. This is an 
effect of the coupling of the Cooperon to the density, 
which is a massless mode at exactly half filling {6 — 0) 
for q = Q Notice that this behavior is peculiar 

of the case p = 1 for a band dispersion arising from a 
nearest-neighbor hopping. In the presence of an attrac- 
tive on-site interaction the system exhibits an enlarged 
symmetry with respect to the instability in the particle- 
particle channel (at q = 0) or in the particle-hole chan- 
nel (at q = Q). Once the symmetry has been explicitly 
broken in the Cooper channel, the density becomes a 
Goldstonc mode, similarly to the phase, reflecting such 
a degeneracy. Since at strong coupling Sp tends to fluc- 
tuate coherently with S\A\ ]l^ , p^ , at small doping and 
strong coupling both XA(q) and Xp{q) have a maximum 
at the wave vector which controls the instability of the 
density mode approaching half filling, i.e., q = Q. 
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FIG. 3. Doping dependence of and given by Eqs. 
and (||) respectively, for t = 1 and U = 20. 




As a result, for < the long-distance behav- 
ior of ^A(r) should be controlled by the characteris- 
tic length obtained by considering the expansion of 
Xa{(i) around Q, 



^A(r) 



dq 



M2+A(q-Q)2 

g^qr 



M2 + Aq2 ' 



(31) 



so that the resulting Xa{j^) is a staggered function with 
an exponential envelope controlled by the stiffness A of 
the q- modes near Q and by the mass . The parameter 
A is obtained by evaluating, at q = Q, the second-order 
derivative of A'^^(q), as given by the expression (p^, 



A^ \dlX^'(ci] 



q=Q 



1 - 26^ 
2U 52 



(32) 



According to Eq. ( pi| ) , the long-distance decay of Xa (r) 
is exponential, with a characteristic length scale 



A 



\ Ml ~ 



1-252 
8^2 ' 



(33) 



which matches continuously at 5 = 1/a/2 with the ^0 
previously defined in Eq. (P9[). 



t=l, U=30 Diluted Limit 



t=l, U=30 Near Half Filling 
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FIG. 4. Absolute value of X/^{r) normalized to its r = 
value for r = nx at various 5. Left panel: the slope of the 
long-distance (exponential) decay of |XA(r]| in the diluted 
limit is given by — CcT^i with from Eq. (E9|). Right panel: 
|^A(r)| near half filling. Although the characteristic length 
of the exponential tail increases as the doping decreases, 
|JfA(r)| is strongly suppressed at much shorter distances. 
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According to our strong-coupling 
^0 and vanish at (5 = 1 / v^. This 
of retaining only the order O(a^) in 
X^^iq). By including in Eqs. @- 
which are of the form ^ a^{w^(q) -|- 
sin''(gy/2)]}, the ratio between the 
the mass term in Eq. gives a 

order Co --C^-Ola'). 



expression ( p4| ) both 
is clearly an artifact 
the q dependence of 
(pO|) the next terms, 
const X [sm'^{qx/2) + 
coefficient of q** and 
coherence length of 
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FIG. 5. X^^(q) at f = 1, (7 = 30 for 5 = 0.005 < ^o? (solid 
line) and (5 = (dashed line), q varies along the diagonal of 
the Brillouin zone. The width of the minimum at Q decreases 
as the doping decreases, vanishing at half filling. Inset: detail 
of the behavior of X^^(q) away from Q. 

One usually defines the correlation length as 

InXA(r) 



Co 



lim 

r — >oo 



According to this definition, the identification between 
^0 the and for (5 < 1/ \/2 would imply a divergent 
at half filling. Such a divergence is surprising, since in 
this regime one would naively expect that the correlation 
length for the amplitude fluctuations attains a value of 
the order of the lattice constant. Indeed this problem is 
only apparent, as we can see by evaluating numerically 
XA(r) near half filling from Eqs. (|l|)-(0). In Fig. | 
we report Xt\ (r) normalized to its r = value at various 
doping, both in the diluted (left panel) and nearly half- 
filled (right panel) limit. In particular, near half filling, 
it can be seen that, although the long-distance behav- 
ior of ^A(r) is controlled by the large length scale ^^r, 
Eq. (^), which effectively diverges as 5 ^ 0, ^A(r) 
is severely suppressed with respect to its r = value at 
much shorter distances, of the order of the lattice spacing. 
A simple criterion to estimate the coherence length con- 
sists in determining the distance at which Xi^(y) has 
reached a fixed percent k of its r = value, say k = O.I, 
m Fig. H One can identify as the minimum value 



between ^„ and the lattice spacing, which is the lower 
bound for in a lattice model. From Fig. ^ it follows 
that the maximum value of at [//< = 30 is about 5 lat- 
tice spacings, and is reached at J = 0.05. At lower doping 
the decay of (r) is much more rapid, even though the 
long-distance tail extends over a greater distance. On 
the other hand, in the low-density regime the coherence 
length is really diverging. In the left panel of Fig. ^ we 
report |XA(r)|/|XA(0)| at high doping. X/^ij) is expo- 
nentially decreasing in this regime, with the characteris- 
tic length given by Eq. (p9|). Indeed, as b increases, the 
slope of the curve, i.e. — CcT^, decreases, and the overall 
decay of ^A(r) extends up to larger distances. 

The different role played by ^o, Eq- dH), and ^jr, Eq. 
( ^ ) , in determining the decay of Xa (r) can be better un- 
derstood by analyzing in greater detail the q dependence 
of ArA(q) approaching half filling. X^^{q) is reported in 
Fig. Q Let us consider the curve for S — 0.005. Even 
though the minimum of A^a^ is at q = Q, the width of 
this minimmn decreases with 5. This refiects the fact 
that, according to Eq. ( p^ and (30), the divergence of 
^TT, Eq. ( |33| ) at half filling arises from an increasing stiff- 
ness A of the q modes near Q with a fixed mass. As 
a consequence near half filling A"A(q) is almost constant 
(and equal to its q = value) except in a narrow region 
near Q, which then contributes to the integral (|3|) only 
at very large distances and with a small weight. 

In the inset of Fig. || it is shown that ai 6 — 0.005 a lo- 
cal minimum of X^^(q) at q = exists besides the global 
minimum at q = Q. Indeed, the second-order derivative 
of X^^ at q = 0, given by C of Eq. ( |28| ) is again positive 
at ^ < 4a^. At exactly 6 — only the q = minimum 
survives, since, according to Eqs. ( p^ ) and (p^, XpA — 
so that the amplitude fluctuations decouple from the den- 
sity fluctuations and X^iq) reduces to the form As 
a consequence, the maximum of A'A(q) is again at zero 
momentum and the coherence length is 



^oiS = 0) = 



CA 



This small value matches continuously with the charac- 
teristic length scale which controls the short-distance de- 
cay of Xa (r) near half filling, and can be extracted form 
the data reported in the right panel of Fig. ^. 

C. The density-density correlation function 

In the last part of this section, to gain some more 
insight into the physics of the interplay between the 
particle-particle and particle-hole channels at strong cou- 
pling, we compare the above results for the pairing cor- 
relation function with the outcomes of the same calcu- 
lations for the density-density correlation function. In 
analogy with the case of XA^q), starting from the same 
Eq. (Ill), we find that the density-density correlation 
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function, defined as Xp(q, fim) = 2 < SpgSp^q >, is 
given by 



1 ^ t ^ 



XpA(g) 

^ - 5^(9) 



(34) 



i.e. the RPA resummation for Xp in the particle-particle 
channel. Evaluating Eq. ([34| ) in the strong-coupling 
regime by means of Eqs. (^8^(|ol) and Eqs. (|2|)-(|3l), 
we find that at leading order in a X~^{q) = M^+Cw{q), 
where M"^ = X-^{0) = Aa'^/U, and C = ia^/2U){2S^ - 
1)/(1 - ^2). It follows that C > at (5 > 1/V2, so that 
the maximum of Xp{q) is at q = and the coherence 
length ^p, which controls the exponential decay of Xp{r), 
does coincide with ^o, Eq. ([29|), 



26^-1 

S(l-J2)- 



(35) 



Instead, for (5 < 1/ \/2, C < and the maximum of Xp{q) 
is at q = Q. Indeed, the mass term for the density- 
density correlation function at q = Q is 



Ml^x;;\Q)=M' 



8C 



4a2 



U 1-52' 

and the second-order derivative at q = Q is 
1 



(36) 



-C. 



q=Q 




As a consequence, for 5 < l/\/2, AlJ < -M^, and the 
function Xp{r) oscillates with an exponential envelope 
controlled by the characteristic length 



(37) 



which coincides with ^^r, Eq- (|33|). However, contrary to 
the case of the |A| mode, which is suppressed at much 
shorter distances due to the presence of a second small 
length scale associated with the secondary minimum of 
X^^(q) at q = 0, the exponential decay of |Xp(r)| is con- 
trolled by a single length scale given by Eq. (p7|), which 
diverges at S = 0. This divergency is due to the fact that 
the density mode becomes a Goldstone mode for the en- 
larged symmetry of the Hubbard model at half filling, so 
that the mass M^, Eq. (|3^), vanishes, while the stiff- 
ness A is finite. As a consequence, while for pairing 
fluctuations is ultimately cut off by the lattice constant, 
the divergence of ^p as according to Eq. (5^ is the 
signal of the degenerate instability at p = 1 |[14|. On the 
other hand, since in the diluted limit the pairing and den- 
sity fluctuations become proportional, as shown in Ref. 
p3[ , the characteristic length of the spatial decay of both 
pairing [Eq. (p9|)] and density mode [Eq. (p^)] diverges 



according to the mapping with the bosonic model. We 
like to point out that the degenerate instability at p = 1 
is specific of the Hubbard model with nearest-neighbor 
hopping. Indeed, as discussed in Appendix B, the inclu- 
sion of a next-to-nearest-neighbor hopping term t' spoils 
this degeneracy and the charge sector becomes massive 
even at p = 1. However, by approaching p — I, the 
charge mode at q = Q still affects the pairing correlation 
function as discussed in this section, provided t'/t is not 
too large. 

The extension of the above results to the strong- 
coupling regime for the d-wave symmetry is not straight- 
forward. Indeed, the decoupling of the interaction Hj, 
Eq. (^), in the particle-hole channel does not simply 
reduce to a density-density interaction, but rather renor- 
malizes the bare band dispersion, limiting the possibility 
for an analytical treatment. Nevertheless, supported also 
by the analysis of the weak-coupling case, we expect that 
in the strong-coupling limit results similar to those of the 
s-wave superconductor hold. In particular, one would 
find a coherence length of the order of the lattice pa- 
rameter for all doping except in the low-density regime, 
where the mapping to the bosonic system should be re- 
covered. Anyway, further numerical work is required to 
address this issue in more detail. 



IV. CONCLUSIONS 

In summary, we studied the evolution of the super- 
conducting coherence length from weak to strong cou- 
pling. For both s-wave and d-wave superconductors, is 
defined through the spatial decay of the correlation func- 
tion for the fluctuations of the modulus of the order 
parameter. At weak coupling the issue arises of the com- 
parison between and the Cooper-pair size ^pair ■ These 
two length scales are both of order vp/^o in s-wave su- 
perconductors. However, as we showed, this identifica- 
tion of ^0 and S,pair at weak coupling does not hold for 
d-wave superconductor, where S^pair is divergent, while 
is finite and of order vf/Aq. 

At strong coupling the modulus and density fluctua- 
tions are coupled, leading to a contribution of the density 
mode to the correlation function Xj\ . The case of s-wave 
pairing within the negative-?/ Hubbard model is simpler 
and allows for a detailed analytical treatment. We then 
evaluated i^o for s-wave superconductors by properly in- 
cluding density fluctuations. We found that the coher- 
ence length in the diluted regime diverges as ~ 
according to the mapping of the Hubbard model into a 
weakly-interacting effective bosonic model, in the strong- 
coupling and low-density regime. Away from the diluted 
limit the coherence length attains a value of order of few 
lattice spacings, as reflecting the local character of the 
Cooper pairing. Similar results are expected in the case 
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of d-wave pairing, where however the analytical treat- 
ment cannot be carried out, and a deeper insight can 
only be gained by further numerical investigation. 

Finally, we comment on the issue of choosing an ap- 
propriate cutoff for the phase-only action when analyzing 
phase-fluctuation effects in superconductors. Our results 
indicate that at weak coupling it is reasonable to consider 
momenta up to a cutoff |qc| — Aq jvp for both s-wave and 
d-wave superconductors. In the strong-coupling regime, 
as we explicitly showed in the s-wave case, one should 
take Iqd — 1/a over a wide doping region, with the no- 
ticeable exception of the low-density region, where di- 
verges and Iqd — ^/p/a. 

Acknowledgments: We thank M. Capone, C. Di Cas- 
tro, M. Grilli, P. Fieri, and G. Strinati for stimulating 
discussions and suggestions. 



APPENDIX A: COOPER-PAIR SIZE AND 
COHERENCE LENGTH AT WEAK COUPLING 



In this appendix we report the detailed calculation of 
(,pair and ^0 in the s-wave and d-wave case at weak cou- 
pling. Let us start with ^pair, which is defined in Eq. (|l|) 



as £,pair = yJ^T/V where 



7V= / dr|^(r 



* k 



(^kAk) • (^k^) , 



(Al) 



and 



k 



In the s-wave case only the second term in the right-hand 
side of Eq. ( |Al[ ) survives, and letting = (9kCk)^ we 
have 



1 , A2;;2 2 
J_ \ " ^0?k^k 

^ s 

An / dx 



An / dx 



(x2 + A2)3 Ns 



V{x), 



where we define V{x) = (1/-/Vs) X^k '^k^(^ ~ Ck)- Since 
at weak coupling the main contribution to the above 
integral comes from a; ~ 0, we take V{x) ~ V{Q) = 



(l/Ns) X^k '^k'^(Ck) and we extend the integral over x be- 
tween — oo and -l-oo to extract the leading behavior 



Similarly, 



dx 



x"^ _ irViO) 
N{x) 



^iV(0)Ao, 



where N{x) = (l/iVs) J2k H^^^^) is the density of states 
in the metallic phase. Then we find 



8AI 



(A2) 



where, according to Eq. (|T^, we introduce the mean 
value of the velocity at the Fermi surface 



N{oy 



It is wor th noting that the approximations leading to 
Eq. (A2) arc no longer valid when the chemical poten- 
tial approaches the extrema or the saddle points of the 
spectrum (critical points), where a more refined eval- 
uation is needed to get the correct result. However, by 
means of numerical calculations we checked that the es- 
timate ( [A^ ) is very accurate even close to the critical 
points. 

In the d-wave case £,pair is infinite, since the nodal 
quasiparticles giv e a l ogarithmically divergent contribu- 
tion to J\f in Eq. ( |Al[ ) . Near the nodes the quasiparticle 
spectrum E^ is cone-like, ~ wifci, Ak — ^2^2, and 
i^k — \/'v\k\ -I- Here vi and V2 are the Fermi ve- 

locity and the slope of the gap Ak at the node, and fci, k2 
are the components of k along the directions perpendic- 
ular and parallel to the Fermi surface respectively, mea- 
sured from the node. By introducing the polar coordi- 
nates {E, 9) , such that when k is near a node ^k = cos 9 
and Ak = -Esin^, and using the identity, 



-T 

fci fc2 



d9dEE 
V1V2 



we find that the first two (positive) terms in the right- 
hand side of Eq. ( Al) are divergent, while the last term, 
which in principle has not a definite sign, vanishes due 
to the orthogonality of 9kAk and 9kCk for a cone-like 
spectrum. For example, the first term gives 



"Hi 

V2 



-y 



k 

2tt 



^(SkAk)^ 



de 



^ ,^E^co8-^9 

dh cx 

£6 



dE 



where A is a properly defined upper cutoff. Notice, in- 
stead, that the contribution of the nodal quasiparticles 
to V is finite. 
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1 



A2 

. El 



APPENDIX B: THE EFFECT OF A 
NEXT-TO-NEAREST-NEIGHBOR HOPPING 



1 



V1V2 Jo 



27T 



d0 



dE 



E^ 



sm 



^2 



dEE. 



We next turn to the evaluation of the coherence length 
^0 in the weak-coupling regime, defined in Eq. (|l^. It 
can be seen that at weak coupling the leading contribu- 
tion to c conies from the last term in Eq. so in Eq. 
(p^ we must take 



5 v^7^ 
16N, ^ El 



while is given by Eq. (||). Observe that both c and 
are finite whatever is the gap symmetry. In the d- 
wave case it can be seen that, thanks to the presence 
of the 7^ factor in the previous equation, the contribu- 
tion of nodal quasiparticles to c is indeed proportional 
to / dOdEE'^ sin* 6 cos^ 6 / E'^ . Let us more generally con- 
sider the case of arbitrary gap modulation 7k and band 
dispersion ^k- We first evaluate the mass term 



m 



A2 
2Ns 



E 

k 



A2 
^0 



dz 



G{z) 



where G(z) = (1/iV.) Ek l7k|<5 (z - Ck/7k). Again, the 
main contribution to the above integral comes from z « 
0, and we obtain 



m 



G(0) 



1 



Analogously, we have 



16 



dz 



iz^ 



A2 
^0 



-W{z) 



W{0) 
12Al' 



where W{z) = (1/iV,) ^^K/lTkO-^ (z - $u/7k) and 
W{0) = V{0). As a consequence, according to the defi- 
nition (\L2\}, we can write 



ViO) 



12A2G(0) 



12 A|' 



(A3) 



which represent the generalization of the standard result 
for a continuum model, to the lattice case with arbitrary 
band dispersion ^k, and for arbitrary symmetry 7k of the 
gap parameter. 

In the s-wave case Eq. (A3) reduces to = u|n/12AQ, 
so that S,pain Eq- (A2), and ^0 differ only by the numer- 
ical factor, ^pair/Co — 1/3/2. Both lengths are, in turn, 
proportional to the Pippard length scale ^e.m. = vp / ttAq, 
which appears in the electromagnetic response function 



Although the strong-couphng results obtained in Sec. 
Ill are specific of the negative-C/ Hubbard model with 
nearest neighbor hopping, we argue that they should not 
be dramatically modified by extensions of the original 
model. In particular, we discuss here in some detail the 
effect of a next-to- nearest- neighbor hopping term t' . It is 
known that such a term breaks the extended symmetry 
of the half-filled model with t' = 0, making charge order- 
ing unfavorable. At weak coupling this is the case, due to 
the spoiling of the perfect nesting of the Fermi surface as 
soon as t' ^ 0. Therefore, at weak coupling the supercon- 
ducting state is favored at all doping down to (5 0, and 
never becomes degenerate to the charge-ordered state. 
Since at small U the coupling of the particle-particle and 
particle-hole channel is negligible, the weak-coupling re- 
sults discussed in Sec. II are unaffected by the modifica- 
tion of the charge sector due to the t' term. 

At strong coupling the extended symmetry is spoiled 
as well by the inclusion of t' and again superconductivity 
is favored with respect to charge ordering, even at p = 1. 
This can be more easily understood exploiting the spin 
analogue of the extended symmetry, which is achieved by 
mapping the negative-C/ Hubbard model onto a positive- 
U model at half filling. We point out that, contrary to 
the case t' = 0, in which the repulsive Hubbard model is 
recovered [ pO[ , the transformed t' term acquires an extra 
sign which depends on the spin of the hopping fermion, 
t'{~U) — > at'{U). Subsequently, in the limit of large 
U, we map this positive-C/ model onto an effective spin 
model which is the usual antiferromagnetic Heisenberg 
model with coupling J ^ for nearest-neighbor 

spins, but yields the non SU (2)-invariant coupling 

J' [s^s]~is:[s^ + sfsy)], 

with J' ~ f^/lUl, between spins located on next-to- 
nearest- neighbor sites << i,j >>. Evidently the mag- 
netic order along the 2;-axis (which corresponds to charge 
ordering in the negative-f/ Hubbard model) is frustrated 
due to the competition of the two antiferromagnetic 
couplings, whereas the magnetic order on the xy-plane 
(which corresponds to superconductivity) is not frus- 
trated due to the absence of interference between an 
antiferromagnetic nearest-neighbor and a ferromagnetic 
next-to-nearest-neighbor coupling. Indeed, the charge 
sector acquires a mass of order t''^/U at half filling. Nev- 
ertheless the results discussed in Sec. HI are not dramat- 
ically modified, as long as the mass in the charge sector 
stays small near half filling, so that the strong peak at 
q = Q survives in Xp{q). 
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